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A mathematical model is developed for the prediction of scattered broadband shock- 
associated noise. Model argnments are dependent on the vector Green’s function of the 
linearized Euler equations, steady Reynolds-averaged Navier-Stokes solutions, and the two- 
point cross-correlation of the equivalent source. The equivalent source is dependent on 
steady Reynolds-averaged Navier-Stokes solutions of the jet flow, that capture the nozzle 
geometry and airframe surface. Contours of the time-averaged streamwise velocity compo- 
nent and turbulent kinetic energy are examined with varying airframe position relative to 
the nozzle exit. Propagation effects are incorporated by approximating the vector Green’s 
function of the linearized Euler equations. This approximation involves the use of ray the- 
ory and an assumption that broadband shock-associated noise is relatively unaffected by 
the refraction of the jet shear layer. A non-dimensional parameter is proposed that quanti- 
fies the changes of the broadband shock-associated noise source with varying jet operating 
condition and airframe position. Scattered broadband shock-associated noise possesses a 
second set of broadband lobes that are due to the effect of scattering. Presented predictions 
demonstrate relatively good agreement compared to a wide variety of measurements. 
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Description 
Prefactor constants 
Anisotropic turbulence coefficients 
Speed of sound 

Coefficients of turbulent scales 

Nozzle exit diameter 

Fully expanded diameter 

Diffraction coefficient 

Vector from source to airframe edge 

Frequency 

Green’s function 

Turbulent kinetic energy or wavenumber 

Axial wavenumber of ps 

Component of turbulent length scale 

Design Mach number 

Fully expanded Mach number 

Pressure 

Shock pressure 

Distance to observer from nozzle exit 
Vector from source to observer 
Vector from airframe edge to observer 
Spectral density of acoustic pressure 
Strouhal number 
Time 

Velocity vector 
Turbulent velocity scale 
Velocity component of vector 
Green’s function of LEE 
Observer vector 
Axial distance from nozzle 
exit to trailing edge 
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Source vector 

Radial distance from nozzle 
exit to trailing edge 
Vector from source to observer 
Wedge angle 

Non-dimensional parameter 

Ratio of specific heats 

Dirac delta function 

Jet spreading rate 

Dissipation 

7t//3 

Pressure component of vector 

Green’s function of LEE 

Density 

Retarded time 

Turbulent time scale 

Polar angle to diffracted ray 

Polar angle relative to incident ray 

Angle to observer from 

upstream nozzle axis 

Frequency 


Abbreviations 

BBSAN Broadband shock-associated noise 

CFD Computational fluid dynamics 

FUN3D Fully-unstructured Navier-Stokes 

LEE Linearized Euler equations 

RANS Reynolds- averaged Navier-Stokes 

SPL Sound pressure level 

TTR Total temperature ratio 
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Introduction 


The fluid dynamics of heated, convecting, compressible turbulence that contains shock waves and Prandtl- 
Meyer expansion waves is not well understood. Interactions between shock waves and turbulence creates 
shock-associated noise. One such canonical flow that produces shock-associated noise is the off-design su- 
personic jet. Off-design supersonic jets are observed in aerospace applications and include the exhaust of 
rockets and air-breathing engines (e.g. jet plumes). In practice, jet engines are integrated with the flight 
vehicle airframe, and this integration can impact the flow and noise. Broadband shock-associated noise (BB- 
SAN), which is one component of jet noise, is scattered by the flight vehicle airframe. This paper presents a 
predominantly first principles approach, relative to others, using an acoustic analogy to predict the scattered 
BBSAN. 

In this paper, we use the partially comprehensive acoustic analogy of Miller^ to predict the BBSAN 
component of jet noise in the presence of an airframe scattering surface. Model arguments are dependent 
on the vector Green’s function of the linearized Euler equations (LEE), steady Reynolds-averaged Navier- 
Stokes (RANS) solutions, and the two-point cross-correlation of the equivalent source. The equivalent 
source is dependent on steady RANS solutions of the jet flow that capture the nozzle geometry and airframe 
surface. The jet mean flow and turbulent statistics are obtained by solving the steady RANS equations 
using a computational fluid dynamics (CFD) approach. Contours of the time-averaged streamwise velocity 
component and turbulent kinetic energy are examined with varying airframe position relative to the nozzle 
exit. Propagation effects are incorporated by approximating the vector Green’s function of the LEE. This 
approximation involves the use of ray theory and an assumption that BBSAN is relatively unaffected by 
the refraction of the jet shear layer. The vector Green’s function of the LEE is written as a function of the 
Green’s function of the Helmholtz equation. The latter Green’s function is calculated with the ray-tracing 
approach of Miller.^ 

This paper first surveys relevant measurements and prediction models for BBSAN, though no review of 
a prediction model for scattered BBSAN is presented, as this model is the first that includes scattering. The 
mathematical model, its associated arguments, and its implementation are described. Emphasis is placed on 
the physical meaning of the model arguments. Details of the approach to And the steady RANS solution are 
presented. A non-dimensional parameter is described to illustrate the deformation of the shock cells within 
the jet plume. A subset of the steady RANS solutions and BBSAN predictions are presented; the latter are 
compared with measurement. Finally, conclusions are drawn from the present theoretical investigation. 

A Survey of Previous Measurements 

A large number of investigators conducted careful measurements to study BBSAN. The first successful 
characterization of BBSAN was likely performed by Harper-Bourne and Fisher.^ They showed that BBSAN 
is caused by the interactions of coherent turbulence with the semi-periodic shock cell structure within the jet 
plume. BBSAN is dominant in intensity relative to other sources at mid- to high-frequencies in the upstream 
and sideline direction relative to the jet flow direction. The spectral content consists of multiple broad peaks 
that decrease in intensity with increasing frequency. Only a few years before the work of Harper-Bourne 
and Fisher,^ Yu^ conducted near-held narrowband acoustic measurements from a wide range of jets. These 
measurements captured the near-field radiating nature of the noise intensity from off-design supersonic jets. 
Similar to the measurements of Yu,^ Tanna et al.® conducted measurements using a convergent nozzle with 
0.0508 m exit diameter and presented a large range of narrowband auto-spectra. The studies of Yu^ and 
Tanna et al.® represented large scale measurement programs that made available large databases of acoustic 
measurements from off-design jets. 

After the work of Yu^ and Tanna et ah,® the National Aeronautics and Space Administration, Langley 
Research Center, conducted a large number of measurements to study BBSAN. In particular. Seiner and 
Norum® A conducted measurements to examine the shock-cell spacing with varying off-design parameter 
and correlated near-field microphone measurements with an intrusive hot-film probe. Comparisons were 
performed between on-design supersonic jets and equivalent off-design shock containing jets. The variation 
of overall sound pressure level (OASPL) between the on-design and off-design cases were quantified. Seiner 
and Yu® conducted near-field correlations between microphones and also have shown correlations between 
early empirical prediction models. Their findings illustrated that the dominant noise source is located near 
the end of the potential core for the conditions they examined. Later, Seiner et al.® examined BBSAN from 
rectangular nozzles. Norum and Seiner^®di conducted measurements of static pressure, and near-field and 
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far-field BBSAN auto-spectra, from unheated jets. The broadband component of shock-associated noise was 
isolated from the discrete component with the novel use of tabs on the nozzle lip. 

Considerable effort was placed on quantifying the effect of forward flight on the BBSAN intensity and 
was an important component for system noise prediction methods. At NASA Langley, Norum and Shearin^^ 
measured the variation of BBSAN in the far-held from under-expanded jets in the presence of hight vehicle 
Mach numbers. Moo < 0.40. They showed that measured BBSAN intensity varies little with varying Moo, 
and prediction techniques such as that devised by Harper-Bourne and Fisher^ were adequate. At higher 
flight vehicle Mach numbers, measurements of Ahuja et al.^^ showed that there is more signihcant BBSAN 
intensity variation. Finally, Yamamoto et al.,^^ in their NASA contractor report, focused on producing 
shadowgraphs from six different nozzle conhgurations with varying jet operating conditions and Moo. 

A number of contemporary measurements have been conducted. These measurements were mainly fo- 
cused on the scaling of BBSAN intensity with stagnation temperature and on the separation of BBSAN 
from the total noise. Viswanathan et al.^^ measured the saturation of BBSAN with stagnation temperature 
and separated the turbulent mixing noise from the BBSAN using incoherent spectral subtraction. Kuo et 
al.^® performed experimental studies to assess scaling laws of BBSAN intensity in the far-held with varying 
jet stagnation pressure and stagnation temperature. Similar to the earlier measurements of Seiner and Yu,® 
Veltin et al.^^ correlated various how-held quantities with the acoustic held but used optical dehectometry. 
Panda^® showed that ‘shock splitting’ occurs within the jet, and is correlated with large organized vortices 
convecting over the shock train. Finally, similar to Norum and Seiner, Andre et al.^® investigated the 
effect of the discrete shock-associated tone on the broadband component. 

Previous Prediction Approaches and their Relation to the Present 

Many researchers created theoretical models for the prediction of BBSAN. Perhaps the hrst successful pre- 
diction method for BBSAN was created by Harper-Bourne and Fisher.® The model was based on the premise 
of coherent interaction between turbulence in the jet shear layer and the nearly periodic jet shock cell struc- 
ture. This interaction was modeled as a series of correlated point sources that combine constructively and 
destructively. Harper-Bourne and Fisher’s® model depended on the rate of decay of turbulence between 
shocks and a characteristic spectrum produced by each shock wave shear layer interaction. These quantities 
were obtained with least square fits to match the far-held data. 

Tam^® developed a stochastic model for BBSAN where the basic physical model was described by Tam 
and Tanna.^® The shock cell structure was modeled as a waveguide following the method of Pack^^ and 
the effect of divergence and turbulent dissipation was included following the method of Tam et al.^® The 
large-scale turbulence in the jet shear layer was modeled as a random superposition of instability waves 
supported by the jet mean flow, as described by Tam and Chen.^^ Tam^® modihed the model of Tam^® 
to predict BBSAN from heated jets with a moderate off-design parameter. Tam and Reddy^® developed a 
method for the prediction of BBSAN from rectangular nozzles with aspect ratios up to six. 

Morris and Miller^^ developed a prediction method for BBSAN that is based upon an acoustic analogy 
approach. The Euler equations were rearranged into the LEE operator and were equated to equivalent 
sources. The convolution of the equivalent sources with the vector Green’s function of the LEE resulted in 
a closed-form analytic expression for the spectral density of the pressure. Equivalent sources were modeled 
based upon a two-point cross-correlation and steady RANS solution. The two-point cross-correlation was 
modeled using dimensional arguments involving the shock pressure, turbulent velocity fluctuations, and 
local scales of turbulence. This formulation was based upon a separable model approach that is similar to 
the model of Ribner.^® The resultant model consisted of a volume integral containing the jet plume and 
an integral of the shock cell pressure wavenumber spectrum. Their acoustic analogy successfully predicted 
BBSAN for a wide range of single-stream jet operating conditions. 

The model of Morris and Miller^^ was extended by Miller and Morris^® to predict BBSAN from dual- 
stream and rectangular jets. Their model made no assumptions regarding the geometry of the nozzle. The 
canonical BBSAN models of Harper-Bourne and Fisher,® Tam,^® and Miller and Morris^® do not account 
for the refraction of sound through the jet shear layer. Miller and Morris®® studied propagation of BBSAN 
by altering the vector Green’s function of the LEE to contain the Green’s function of Lilley’s equation as 
an argument. Henry et al.®® studied propagation of BBSAN by using a ray tracing approach. The Green’s 
function of Lilley’s equation within Miller and Morris®® was found numerically using an adjoint approach 
and a locally parallel flow assumption. It was shown that in the upstream and sideline radiation direction 
refraction effects have little effect on BBSAN. In the downstream direction refraction has a larger effect but 
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the source of BBSAN is relatively less intense. However, a small improvement in prediction accuracy was 
shown relative to predictions excluding refraction effects. 

As jet stagnation temperature increases the resultant intensity of BBSAN increases then ‘saturates’ 
(ceases to increase with increasing temperature) as observed experimentally by Viswanathan et al.^® The 
model of Harper-Bourne and Fisher^ and Morris and Miller^^ do not account for varying jet stagnation tem- 
perature while the model of Tam^® does. This latter model accounts for slight decreases of BBSAN intensity 
at all angles with increasing temperature using a simple empirical correction factor that is proportional 
to the inverse of the fully expanded density. Miller^^ reformulated the model of Morris and Miller^^ and 
included BBSAN refraction effects using the approach of Miller and Morris^® to account for BBSAN satura- 
tion with increasing temperature. Miller^^ accurately predicted the saturation of BBSAN and proposed that 
saturation is due to an exact balance between increasing source strength and decreasing refraction intensity 
amplification. 

Very recently, Miller^ created an acoustic analogy for jet noise that included equivalent sources for 
both turbulent mixing noise and BBSAN. It was shown to predict the total noise accurately for a wide 
range of jet Mach numbers and jet stagnation temperatures. The turbulent mixing noise and BBSAN were 
also predicted independently. The model made no assumptions regarding fine-scale or large-scale similarity 
spectra or sources, shear-noise or self-noise equivalent sources, or other fictitious sources. 

A scattering approach for jet mixing noise was recently developed by Miller.^ They adopted the acoustic 
analogy of Miller^^ that was created for the prediction of jet noise ground effects. The Green’s function of 
Lilley’s equation was approximated by creating a bridging function of asymptotic solutions. It contained 
arguments involving the locally parallel mean flow and the Green’s function of the convective Helmholtz 
equation. A ray tracing solver, for the Green’s function of the convective wave equation, was used to find 
this latter argument of the bridging function. Accurate predictions of jet mixing noise from jets operating 
over a wide range of Mach numbers and temperature ratios agreed with the measurements of Brown. Like 
the approach created for scattered jet mixing noise of Miller^^ and Miller,^ the prediction method presented 
here is complimentary, but for scattered BBSAN. 


Mathematical Approach 


The predictions shown within this paper are based upon the partially comprehensive acoustic analogy of 
Miller.^ The far- field model equation contains equivalent source terms for the mixing and BBSAN compo- 
nents of jet noise. Details of the developed equation are shown in Miller^ and are heavily based upon the 
theories of Tam and Auriault,^^ Morris and Farassat,^® and Morris and Miller. The model is dependent on 
steady RANS solutions and an analytical model of the two-point cross-correlation of the equivalent sources 
of mixing noise and BBSAN. The spectral density, S, at observer location, x, is 
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where the integrand is a combination of spatially varying equivalent sources of mixing noise, BBSAN, and 
the vector Green’s function of the LEE. The integral is restricted to the jet plume. The observer vector x 
is a vector from the nozzle exit to the observer, and the source vector y is a vector from the nozzle exit to 
the sources. Note that bold symbols denote vectors and a non-bold counterpart represents the magnitude 
of the vector. The angular frequency is w = 27 t/, where / is the frequency. Spatially varying time-averaged 
quantities consist of the velocity, u = m(m, u, tc), and density, p. Statistical quantities related to the RANS 
closure are the turbulent kinetic energy, fc, and specific dissipation, e. Turbulent scales, that are related 
to the RANS solution, are the streamwise, lx-, and cross-stream lengths scales, ly and /^, the velocity Ug, 
and time scale, Tg. The coefficient matrix amn describes anisotropic effects on the BBSAN component of 
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the total noise and can potentially be described following the work of Khavaran.^^ Here, we let amn be 
unity. Ambient reference values are denoted by subscript oo. The shock pressure, Ps = p — Poo, represents 
the pressure due to the shock cell structure. Within the inner integrand, the wavenumber spectrum of Ps is 
denoted by Ps, and ki is the wavenumber. 

Though Eqn. f is relatively complicated, individual terms have physical meaning. The first and second 
terms within the integrand are a result of analytical integration of a model for the two-point cross-correlation 
of the equivalent sources within the jet. For simplicity, it is assumed that this correlation has the same 
functional form for all equivalent sources. The first term within the integrand controls scaling of intensity 
for all sources with respect to the turbulent scales, u, and oj. The second term within the integrand is a direct 
result of the mathematical form of the two-point cross-correlation, that (in this model) is assumed Gaussian 
for simplicity. Though, a better model would likely involve a combination of Gaussian and exponential decay 
with separation distance in temporal, radial, and axial directions, this approach only assumes the former. 
The first and second terms within the braces correspond to the mixing noise generated by fluid dilatation 
and unsteady forces per unit volume. Together they predict the total mixing noise produced by the jet. The 
last term within the braces of the integrand is the source term of BBSAN. BBSAN intensity is proportional 
to the turbulent kinetic energy, k, and the shock pressure, Ps, squared. The integral involving ki reduces to 
the shock pressure within the jet. 

The turbulent scales are related to a steady RANS solution by, lx = cik^^'^e~^, Tg = Crke~^, and 
Us = Cu(2fc/3)^/^, where c;, Cr, and Cu, are constant coefficients. These coefficients are calibrated with 
an over- and under-expanded jet at the sideline observer location. The under-expanded jet operates at 
Mj = 1.50 and total temperature ratio, TTR = 1.00, from a convergent nozzle with exit diameter of 0.0508 
m. The over-expanded jet operates at Mj = 1.29 and TTR = 1 from a Md = 1.50 nozzle with exit diameter 
of 0.0508 m. Once the coefficients are obtained they are held constant regardless of nozzle geometry or 
jet operating condition. Also, we let Cs equal unity. The coefficients are Cr = 1.25 and c; = 1.25 for over- 
expanded flows and Cr = 0.50 and ci = 3.00 for under-expanded flows. The cross-stream length scales are 1/3 
of the streamwise length scale. The convective velocity or Mach number is approximated as 7/10 the local 
value. Frequency dependence of the length scale has shown to improve predictions relative to measurement, 
especially for mixing noise, but demonstrates negligible improvement for BBSAN predictions. It could be 
included by following the approach of Lieb and Goldstein.^® 

Within Eqn. 1, represents the or logarithmic component of pressure. These components are 
governed by the solution of a set of partial differential equations 


and 


dt 


- 


dxi 


6{x-y)S{t-T) 6on 


( 2 ) 


dt dxj 


+ v" 


dvj 

dxi 


cItt" 

dxi 


6{x-y)6{t-T) Sin, 


( 3 ) 


where represents the velocity component. The Dirac delta function is S and Sij = 1 for i = j and 0 for 
i ^ j. Assuming that the flow is statistically stationary, we transform tt^ from the time to frequency domain 
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where r is the retarded time. 

This mathematical survey is focused on BBSAN and not the total noise or mixing noise (see Miller^ 
for details), thus we focus on isolating the BBSAN component of the total noise from Eqn. 1. First, the 
coefficients within Eqn. 1, Ag and Bg, are equated to zero and after simplification Eqn. 1 is 


/ OO pOO p 

/ / 

-OO J — oo J — 


T^^^'^Cjp^d^^TT^klxlylzPgX'^Ts 


/_oo + (UXI + VX2 + WX3 + CaoX)^T^0J^) 

/ OO 

Ps exp[-ikiyi]dkidy , 

-OO 


( 5 ) 


5 of 24 


American Institute of Aeronautics and Astronautics 


where summations of m and n are implied from one to three and arguments of the terms are dropped for 
compactness. Equation 5 is evaluated by approximating the spatial integrals numerically. The majority of 
the arguments are dependent on the steady RANS solution in the vicinity of the airframe over the volume 
of the jet plume. The wavenumber of the shock pressure, Ps, is evaluated by performing a discrete Fourier 
transform in the direction of the flow. Integration of Ps exp[— j/cij/i] is performed numerically (without the 
use of a Fourier transform library), and the limits of integration are restricted from |fei| > c, where c > 0 
and represents a small constant. This approach removes the ‘DC’ component of the wavenumber frequency 
spectrum and has the advantage of retaining multiple dominant peaks of Ps- This approximate approach for 
evaluating the inner integral is beneficial because of the use of approximating the shock pressure, Ps, about 
Poo ■ A superior approach might use a shock- free base flow. 


The Scattered and Quiescent Field 

This paper compares the predicted BBSAN that is scattered by the airframe to the predicted BBSAN that 
propagates freely. This comparison requires that the steady RANS solution be evaluated both with and 
without the presence of the airframe. The presence of the airframe alters the statistics of the steady RANS 
solution, and these terms reside in the integrand of Eqn. 5. This approach captures the airframe’s effect on 
the equivalent source. 

The term, 7r*"7r™, in Eqn. 5 alters the intensity due to propagation effects. Thus, a tailored form of 
7Tg captures the effect of the scattered field. Recall that jet shear layer refraction effects have little impact 
on the variation of far-held BBSAN auto-spectra. We adopt the approach of approximating the meanhow 
as either uniform or quiescent for propagation purposes only. In the absence of a mean how (the source 
terms are still dependent on the steady RANS solution) Eqns. 2 and 3 are simplihed, performing the Fourier 
transform with respect to time, and solving for tt”, results in the vector Green’s function of the LEE for a 
quiescent environment 


T^g{x;y,uj) = -^ f -^5{z-y)g{x,z,u)dz, 

Cqo J z OZi 

where 2 ; is a vector from source locations to observers. For free jet predictions 

exp - 2 :|/coo] 

Substituting Eqn. 7 into Eqn. 6 and simplifying yields an expression for the vector Green’s function in 
the absence of varying mean how (which is very similar to the Green’s function of the Helmholtz equation) 
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When a scattering surface is placed in the held of the jet plume, a tailored form of must be found. We 
approximate the tailored form of by hnding a tailored Green’s function of g. Recall that BBSAN occurs 
at mid- to high-frequencies. Thus, one fast approach to calculate g in the presence of a scattering surface 
is with the use of geometric acoustics. Recently, Miller^ showed that a ray tracing approach can be used to 
accurately predict jet mixing noise for mid- to high-frequencies. We use the same ray tracing approach for 
these predictions, and expect this approach to be a better approximation than previously, due to the fact 
that BBSAN occurs mainly at mid- to high-frequencies. The tailored Green’s function is approximated as 
the sum of the incident and diffracted field 


g{x;y,u}) = gi{x;y,u}) + ^gD{x-,y,uj), (9) 

i=l 

where gj is the incident ray, equivalent to the free-held Green’s function, and gu are the scattered rays. If 
the observer is in a shadow region, the incident ray is zero and the total diffracted held is the sum of all 
diffracted rays. Keller^® observed from Sommerfeld’s^° solution of diffraction by a semi-inhnite screen with a 
straight edge that incident waves propagating in a direction normal to the edge create a cylindrical diffracted 
wave centered on the edge. That is, when acoustic rays dihract they have similar properties to the initial 
incident rays but depend on the point of diffraction. Assuming that the amplitude of the diffracted ray is 
proportional to the amplitude at the point of dihraction, the acoustic pressure of the diffracted ray is 
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p{x,uj) = Dp{E,uj) 
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where D is the diffraction coefficient and ri = x — E. The acoustic pressure field at the diffraction point 
on the edge is, p{E,uj) = Aexp[ikE] \E — y \~^ , where E is the diffraction location on the edge and A is a 
constant. Keller^® showed that Eqn. 10 satisfies Sommerfeld’s^° exact asymptotic high-frequency solution 
for diffraction of a wave by a half-plane. After adopting a simplification from Agarwal et al.^^ the resultant 
diffraction coefficient is 
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where (ps is the polar angle to the incident ray, (j)o is the polar angle to the diffracted ray, 6 is the oblique 
angle between the edge and incident ray, v = tt//?, and j3 is the wedge angle. The contribution to go from a 
single diffracted ray is calculated using Eqns. 10 and 11. All diffracted ray contributions are used to compute 
the tailored Green’s function, g. 


Quantification with a Non-Dimensional Number 

We now turn our attention to the development of a non-dimensional number that quantifies the effect of 
the airframe on the BBSAN source. Parameters are identified that alter the noise source with varying jet 
operating condition and airframe position relative to the nozzle exit. 

We propose a non-dimensional number T 

_ DjXp tan[d^] 

-yp{yp-Dl2Y ^ ^ 

where D is the nozzle diameter, Dj is the fully expanded diameter, Xp and yp are axial and radial vector 
components that form a vector from the nozzle exit to the trailing edge, and 5p is the spreading rate of the jet. 
Physically, Eqn. 12 is the ratio of the product of jet and airframe length scales divided by the cross-stream 
length scale and interaction distance. Small values of P imply that airframe effects on the BBSAN sources 
are negligible. Likewise, large values of P imply that the effects of the airframe on the BBSAN sources are 
very large. Note that yp > D/2. 

The fully expanded diameter is 


Dj=D 


/ 1 + ^A/| \ 


( 7 + 1 )/ 4 ( 7 - 1 ) 



(13) 


where 7 is the ratio of specific heats. Values of are not readily available without numerical calculations or 
measurement. The spreading rate is used to estimate the jet impingement location (used in the development 
of Eqn. 12). The jet impingement location is estimated as 


Vp - D/2 
tan [(5^] 


Additional details regarding the development and significance of P are discussed by Miller.^ 


(14) 


Results 

The prediction method outlined in the previous section is validated with a series of measurements con- 
ducted by Brown. These measurements were performed with the Small Hot Jet Acoustic Rig (SHJAR) at 
NASA using the small metal chevron (SMC) series of nozzles. Large flat plates, described by Brown, were 
placed near various SMC nozzles. A photograph of the experiment is shown in Fig. 1 courtesy of Podboy."'^ 
The nozzle and large fiat plate can be seen on the right side of the photograph. The nozzle is pointed 
towards the camera and the large flat plate is mounted vertically. An array of microphones is shown on 
stands centered about the nozzle exit. Note that the plate is between the nozzle exit and the microphone 
array. 
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The experiment of Brown^^ is summarized in the diagram of Fig. 2. The axial distance from the nozzle 
exit to the trailing edge of the plate is, Xp, and the radial distance from the nozzle centerline to the plate 
is, Up. These two quantities and the jet operating condition were varied by Brown^^ to study the scattered 
jet noise. In the measurements the polar microphone array distance was R/ D = 75, and the results have 
been scaled to R/D = 100 for comparison with predictions. Additional corrections were made to account 
for atmospheric absorption. The observer angle, 'h, is measured from the jet upstream axis. For full details 
of the measurement program see Bridges and Brown, Podboy,^^ and Brown. 

Steady Reynolds- Averaged Navier-Stokes Solutions 

The steady RANS solutions are calculated numerically using the fully-unstructured Navier-Stokes (FUN3D) 
solver developed at NASA Langley Research Center. For more information on FUN3D see Anderson and 
Bonhaus^^ or Nielsen. We focus on the SMC016 nozzle, which is a convergent-divergent nozzle with design 
Mach number, Md = 1.50, and exit diameter, D = 0.0508 m. The computational domain extends 1001? 
downstream from the nozzle exit, 501? in the radial direction from the centerline, and 51? upstream from the 
nozzle inlet. Numerical solutions are governed by the steady RANS equations closed by the Menter shear 
stress transport (SST)^® turbulence model. Roe^^ flux vector construction is used for spatial discretization 
and is second order accurate. A Roe^® ‘modmin’ flux limiter is used in cases that fail to converge due to 
very large fluxes in the initial transient solution. The flow-field is initialized as quiescent at the start of the 
simulation. The validation of the FUN3D steady RANS solutions for single-stream jets was performed in an 
earlier study by Smith and Miller. 

One steady RANS solution is shown in Fig. 3 for a jet operating at Mj = 1.29 and TTR = 1. The nozzle 
is convergent-divergent, designed with the method of characteristics, with design Mach number, Md = 1.50 
and D = 0.0508 m. The computational domain contains the nozzle geometry and a large flat plate located 
at Xp/D = 10 and yp/D = —2. Reference Fig. 2 for the dehnition of Xp and yp. For these jet conditions and 
plate position the non-dimensional number, F = 0.309, signifies that the BBSAN source is deformed relatively 
little. Note that the simulations are fully three-dimensional, and these figures only illustrate a plane through 
the nozzle centerline axis and normal to the plate surface. Figure 3(a) shows the unstructured-structured 
computational domain in the x-y plane. The jet plume has been discretized using a cartesian grid and the 
exterior region has been discretized using an unstructured grid. The x-axis and y-axis have been normalized 
by D. The nozzle contour, with nozzle exit located at the coordinate system origin, and flat plate terminating 
at XpjD = 10 and yp/D = —2, are shown in each sub-hgure. There are ample grid points to resolve the 
shock-cell structure within the jet plume, and the boundary layer is resolved along the plate surface facing 
the jet flow. 

Figures 3(b), 3(c), and 3(d) show contours of p, u, and k, respectively. In this particular case, it is clear 
that the entire shock-cell structure is shielded by the airframe surface relative to observers in the sideline and 
upstream direction. The field- variables, p and u, shown in Figs. 3(b) and 3(c), have little deformation due 
to the airframe surface. Note that the contour lines relative darkness corresponds to the subfigure legend. 
That is, they are nearly symmetric about the x-z plane. This is in contrast with the turbulent kinetic energy 
in Fig. 3(d), that does show significant asymmetry, especially so relative to ~p and u. Here, the relatively 
weak deformation of the acoustic source intensity and location caused by the presence of the airframe are 
captured by the steady RANS solution. 

We now examine the effect of the position of the airframe surface relative to the Mj = 1.29 and TTR 
= 1 jet. Figure 4 shows contours of u for Xp = 4 and yp/D = —6, -4, -2, and -1. Values of F are 0.011, 
0.027, 0.124, and 0.741 for Figs. 4(a), 4(b), 4(c), and 4(d), respectively. Contours within Figs. 4(a) and 4(b) 
are nearly identical, while in Figs. 4(c) and 4(d), the potential core length has slightly increased. Though, 
it appears that the flow is still relatively symmetric about the x-z plane. Miller^ showed that as the jet 
approaches the surface, the entrainment is affected between the jet flow and the airframe, and there is 
additional entrainment acceleration compared to the unshielded side of the jet. 

Figures 5(a), 5(b), 5(c), and 5(d) show contours of k for the same jet condition and airframe position as 
Fig. 4. Values of F are the same. Contours of k are very similar between yp/D = —4 and -6, and change little 
when examining steady RANS solutions calculated farther from the surface. Unlike the steady contours of u, 
significant asymmetry is observed at yp/D = —2 and -1, in Figs. 5(c) and 5(d). The amount of production 
of k between the jet and airframe, for yp/D = —2, is signihcantly decreased while the contours of k on the 
opposite side are very similar to the undisturbed cases. This is due to the decreased relative shear layer 
velocity between the jet and airframe as observed in Fig. 4(c). Figure 5(d) shows an opposite trend relative 
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to Fig. 5(c), as there are additional levels of k between the jet and airframe. The additional levels are due 
to the direct interaction of part of the jet plume with the airframe surface. 

We now examine the steady RANS solutions of u and k for the same conditions but extend the airframe 
surface to XpjD = 10. The results are shown in Fig. 6 for u and Fig. 7 for k. Values of F for yp/D = —6 
through yp/D = —1 are 0.028, 0.066, 0.309, and 1.853 respectively. Variation of u in Fig. 6 is very small 
except for yp/D = —1, as shown in Fig. 6(d), where the jet shows a very weak Coanda like effect. This can 
likely be attributed to the additional plane-form of the plate relative to the Xp/D = —4 case. Variation of 
k with yp/D is shown in Fig. 7, and at yp/D = —4 and -6 little difference is observed. Larger asymmetries 
about the x-z plane are observed in Figs. 7(c) and 7(d). These changes in k are much larger than the 
Xp/D = 4 case shown previously in Fig. 5. The magnitude of k in the region of y/D < 0 is lower than 
y/D > 0. The BBSAN source is modihed as it scales with k, especially so on the shielded side of the jet 
flow. Larger values of F that are approximately unity or greater indicate this modification. 

A final contour set of u and k are shown in Figs. 8 and 9 respectively. The plate is extended to Xp/D = 20 
and the other conditions remain the same. Values of F for yp/D = —6 through yp/D = —1 are 0.056, 0.132, 
0.618, and 3.706 respectively. These predictions are very similar to those presented in Figs. 6 and 7. We 
observe the same Coanda effect (an attachment of the jet flow upon the airframe surface) in Fig. 8(d) as we 
did in Fig. 6(d). The extension of the plate by an additional lOH has caused little discernible difference in u 
for this supersonic jet flow. There are some differences in the contours of k between Xp/D = 10 and 20, but 
the contours follow the same trend. In particular, the peak k in the shear layer between the jet and airframe 
has moved closer to the nozzle exit and decreased in magnitude. Generally, the deformation of contours of 
u are relatively small, even for plate positions of yp/D = — 1. Contours of k are much more sensitive to the 
position of the plate, and the importance of intensity of k is shown in the equivalent source model, where 
the BBSAN intensity is proportional to k. The deformation of k for the cases examined is correlated with F. 

Broadband Shock- Associated Noise Predictions 

Here, we will present predictions of scattered and isolated BBSAN. Figure 10 shows predictions compared 
with measurements of acoustic spectra produced by a Mj = 1.29 and TTR = 1 off-design jet. The design 
Mach number is 1.50 and the exit diameter is 0.0508 m. The y-axis represents sound pressure level (SPL) 
per unit St, where St is the Strouhal number, St = fDjuJ^. The observer angle, \k, is shown in the 
upper right corner of each figure and observers are lOOH from the nozzle exit. The plate edge is located 
at Xp/D = 4 and yp/D = —2. The solid and dashed lines represent the free (Bridges and Brown"'^) and 
scattered (Brown^^) jet measurements, respectively. Lines with circles and triangles represent predictions 
for free turbulent mixing noise and BBSAN, respectively, using the method of Miller.^ Predicted scattered 
jet mixing noise is represented by dashed lines with diamonds and was performed by Miller.^ Predictions 
of jet surface interaction noise of Miller^'^ are shown as a line with squares. These predictions are shown 
to illustrate the different components of the total jet noise for the scattered and isolated jet conhgurations. 
Some of these components have been isolated experimentally using the approach of Viswanathan,^^ and are 
similar to present predictions. A description of the physics of these components relative to the presented 
predictions and measurements is outside the scope of this paper. The focus of this paper is scattered BBSAN 
as predicted by Eqn. 5 that is represented by dashed lines with grads. 

Within Fig. 10, except at 4'° = 50, there is very little change in the measured isolated and scattered 
BBSAN. The predictions show larger amounts of shielding relative to measurement at higher frequencies. 
The relative intensity at the peak frequency of BBSAN is captured accurately, as can be seen at 4>° = 50. 
Also, the peak frequency of the scattered BBSAN is higher than the isolated case in this direction. Here, 
F = 0.124 and the steady RANS flow-field is shown in Figs. 4(c) and 5(c). This relatively small value of F 
signifies that the BBSAN source is relatively unaffected by the airframe, and the relative differences between 
the scattered and isolated BBSAN predictions are mainly due to propagation effects. 

Figures 11 through 14 show predictions and measurements produced by the Mj = 1.29 and TTR = 1 off- 
design jet, but with varying yp/D locations. The values of yp/D are -4, -6, -8, and -10, with corresponding 
values of F = 0.027, 0.011, 0.006, and 0.004. Time-averaged steady RANS solutions for these cases are 
shown in Fig. 4 for u and Fig. 5 for k. As with the yp/D = -2 case, these values are relatively small and 
the variation in BBSAN between the scattered and isolated case is due to propagation. At 4>° = 90 and 
110, no measured difference is observed between the scattered and isolated BBSAN. This is likely due to the 
fact that the dominant sound producing shock wave shear layer interactions are not shielded relative to the 
observer. Recall that the dominant source is generally near the end of the potential core, and this region is 

9 of 24 


American Institute of Aeronautics and Astronautics 


not shielded by the airframe surface. Prediction and measurement for the scattered and isolated BBSAN are 
sometimes identical, especially when the observer is not shielded and yp/D is large, which gives credibility to 
the developed approach. It is observed that BBSAN is increasingly shielded and its peak frequency increases 
while = 50 and Xp/D = 4 and the magnitude of Vp/D increase. Predictions capture this trend observed 
in intensity and frequency. By holding dt® = 70 and Xp/D = 4 constant while increasing yp/D, we initially 
observe that the peak magnitudes of scattered BBSAN and isolated BBSAN are in general agreement in 
prediction and measurement. At yp/D = —6, a slight change occurs, and then for yp/D > —8, a large 
reduction of scattered BBSAN is observed. Predictions capture this ‘delayed’ reduction of scattered BBSAN 
intensity between yp/D = —6 and yp/D < —8. 

For the next case, the same jet operating conditions are retained, Mj = 1.29 and TTR = 1, and the plate 
is extended to Xp/D = 10. Predictions and corresponding measurements are shown in Figs. 15 through 18. 
The values of yp/D are -2, -4, -6, and -10, respectively. Values of F are 0.309, 0.066, 0.028, and 0.010. 
Time-averaged steady BANS solutions for these cases are shown in Fig. 6 for u and Fig. 7 for k. Unlike 
the predictions shown for Xp/D = 4, these predictions and measurements demonstrate significant shock- 
associated noise shielding at all angles examined. This is purely due to the extension of the plate by an 
additional 6D and resulting propagation effects. The jet core length is approximately 8D and the dominant 
BBSAN source oscillates in the region between approximately 6D and 9D. By examining Fig. 6, for example, 
for the observer angles chosen it can be observed that the propagation path is restricted to diffracting waves, 
and all predictions and measurements are shielded. The relatively small values of F imply that the source 
of BBSAN is relatively unaffected by the airframe surface, and at most mildly for yp/D = -2 and F = 0.309. 
Very little difference in relative BBSAN shielding is observed between yp/D = -2 and yp/D = -4, even though 
there is a relatively more significant change in F, as shown in Figs. 15 and 16. The quality of the scattered 
BBSAN predictions vary. Generally, higher frequencies are under-predicted, thus corresponding AdB would 
be over-predicted. Though, this would be less of a problem for full-scale engines, because they can exhibit 
greater BBSAN intensity decay at high frequencies. 

The final set of comparisons for the Mj = 1.29 and TTR = 1 jet are shown for Xp/D = 20, and all 
other conditions are held constant. Predictions and corresponding measurements are presented in Figs. 19 
through 21. The values of yp/D are -2, -6, and -10, respectively. Values of F are 0.618, 0.056, and 0.020 for 
yp/D = -2, -6, and -10. Time-averaged steady RANS solutions for these cases are shown in Fig. 8 for u and 
Fig. 9 for k. Like the predictions and measurements for Xp/ D = 10, these predictions also show significant 
airframe shielding at all observer angles and values of yp/D. The scattered BBSAN predictions are more 
aligned with measurement than previous cases. This is likely due to the ray paths having relatively small 
angles relative to the airframe surface. This increases the accuracy of the ray tracing approach, as shown by 
Miller.^ For this reason, near St ~ 10, predictions are in better agreement relative to previous cases, such 
as Xp/D = 4. Though, in the downstream direction, peak scattered BBSAN frequencies are over-predicted 
relative to measurement by St ~ 2. For the yp/D = -6 and -10 cases, shown in Figs. 20 and 21 respectively, 
values of F are very small, and relative error in the predictions and measurement of isolated and scattered 
BBSAN are purely due to propagation. Relative errors between scattered and isolated BBSAN predictions 
and measurement, in Fig. 19 where yp/D = —2, are due to a combination of the ray tracing calculation and 
a slightly deformed BBSAN source. Though, the deformation of the BBSAN source by the airframe surface 
is captured through the steady RANS solution shown in Figs. 8(c) and 9(c). This analysis is complicated 
by the varying intensity of scattered jet mixing noise with varying yp/D, that accounts for some of the total 
measured intensity, relative to the predicted scattered BBSAN intensity. 

Perhaps one of the most interesting features of scattered BBSAN (that was not apparent in the initial 
development of the theory), are elegant ‘double broad lobes,’ that occur in certain circumstances. Recall 
that isolated BBSAN is characterized by multiple broad lobes, that decrease in intensity and spectral width 
with increasing frequency, that form due to constructive interference. Miller®^ isolated these individual 
contributions to the total isolated BBSAN, in the context of the model of Morris and Miller^^ and Tam.^° 
Scattered BBSAN, varying with frequency, is characterized by additional broad lobes at low frequencies, 
that increase in intensity and spectral width below the peak frequency, then decrease in intensity and 
spectral width above the peak frequency. Some examples of this unique scattered BBSAN spectrum can 
be observed in Fig. 13 (xp/D = 4, yp/D = —8, and 4'° = 70 and 90), Fig. 16 {xp/D = 10, yp/D = —4, 
and T® = 90 and 110), or Fig. 21 {xp/D = 20, yp/D = —10, and 4'° = 110). Spectral lobes below the 
peak scattered BBSAN frequency are due to scattering effects. In a mathematical context, they are due 
to the presence of Fresnel integrals that are implicitly contained in analytical solutions of point sources 
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about semi-infinite fiat plates. Spectral lobes at and above the peak frequency are mainly due to traditional 
BBSAN constructive interference and partially due to scattering effects. Above the peak scattered BBSAN 
frequency, the combination of traditional constructive interference and interference due to scattering can 
have a strengthening or smoothing effect on the broadband noise. Nonetheless, the double broad lobes are 
an interesting phenomenon that is unique to shock-associated noise. 

Conclusion 

A mathematical model is developed for the prediction of scattered BBSAN and is based on the partially 
comprehensive acoustic analogy. Predictions using the model, for the scattered and isolated BBSAN, are 
compared with a wide variety of measurements involving a large fiat plate. These predictions are generally 
in satisfactory agreement with measurement. The model arguments are dependent on the vector Green’s 
function of the LEE, steady RANS solutions, and a two-point cross-correlation of the equivalent source. 

The equivalent source is dependent on steady RANS solutions of the jet flow that contain the effects of 
the nozzle geometry and airframe surface. Steady RANS solutions are calculated for an off-design jet with 
various plate locations. Contours of the time-averaged streamwise velocity component and turbulent kinetic 
energy are examined with varying plate position. Generally, the turbulent statistics are more sensitive to the 
airframe than density, velocity, and pressure components. This is due to the sensitivity of turbulent kinetic 
energy on the shear layer gradient that changes due to the variation of entrainment velocity between the jet 
and airframe. 

Propagation effects are captured by finding an approximation of the vector Green’s function of the LEE. 
This approximation involves the use of ray theory and an assumption that BBSAN is relatively unaffected by 
the refraction of the jet shear layer. Also, the dominant energy in scattered BBSAN spectra are at relatively 
mid- to high-frequencies, and the ray theory assumption is accurate in these frequency bands. 

Predictions of scattered and isolated BBSAN agree with measurements of equivalent scattered and iso- 
lated BBSAN. Generally, scattered BBSAN exhibits slightly increased peak frequencies. Also, shielded 
BBSAN exhibits greatly decreased intensities, and generally the model captures these peak intensity and 
peak frequency trends. The model under-predicts very high frequency scattered BBSAN, but the additional 
intensity from scattered turbulent mixing noise will increase the total intensity. 

A non-dimensional parameter, P, is presented that quantifies the changes of the BBSAN source with 
varying jet operating conditions and airframe position. The steady RANS solutions, in particular the vari- 
ation of turbulent kinetic energy with plate position, demonstrate a correlation between P and turbulent 
kinetic energy. Acoustic predictions corresponding to a wide range of P, demonstrate that the model can 
predict the scattering of BBSAN, even with sources altered by the airframe. That is, the model predicts 
scattered BBSAN when the source changes due to the presence of an airframe. It is more likely that the 
approximation of the vector Green’s function of the LEE with the use of ray theory causes more error in 
prediction than the choice of the equivalent source of BBSAN. 

Most interestingly, double broad lobes are observed in the predicted scattered BBSAN. Below the peak 
frequency of BBSAN the lobes are due to diffraction, and above the peak frequency the lobes are due to 
traditional constructive interference. The lobes decrease in intensity as the norm of frequency increases 
relative to the peak frequency. It is likely that the traditional constructive interference, that occurs at high 
frequencies, is disturbed by the constructive contribution due to diffraction. 
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Figure 1. A photograph of the jet structure interaction test courtesy of Podboy.' 










X (Observer) 


Figure 2. Coordinate system of the jet structure interaction test. 





(d) Contours of turbulent kinetic energy. 


Figure 3. Quantities of the steady R.ANS solution on the x-y plane corresponding to the Mj = 1.29 and TTR = 1 jet 
from the S1VIC016 {Md — 1.50) nozzle with D — 0.0508 m. The plate edge is located at XpfD — 10 and ypjD — —2. 
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Figure 4. The Md — 1.50 nozzle operates at Mj — 1.29 and TTR = 1.00. Contours of the time-averaged streamwise 
velocity component are shown on the z — 0 plane. The plate is positioned at XpjD — 4 and at various radial positions, 

Vp!D. 




(c) The plate positioned at ypjD = —2. 


(d) The plate positioned at ypjD = — 1. 


Figure 5. The Md — 1.50 nozzle operates at Mj — 1.29 and TTR, = 1.00. Contours of turbulent kinetic energy are 
shown on the z — 0 plane. The plate is positioned at XpjD — 4 and at various radial positions, ypjD. 
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Figure 6. The Md — 1.50 nozzle operates at Mj — 1.29 and TTR = 1.00. Contours of the time-averaged streamwise 
velocity component are shown on the z — 0 plane. The plate is positioned at XpjD — 10 and at various radial positions, 

Vp!D. 




(a) The plate positioned at ypjD = —6. 


(b) The plate positioned at ypID = —4. 



Figure 7. The Md — 1.50 nozzle operates at Mj — 1.29 and TTR, = 1.00. Contours of turbulent kinetic energy are 
shown on the z — 0 plane. The plate is positioned at XpjD — 10 and at various radial positions, ypjD. 
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Figure 8. The Md — 1.50 nozzle operates at Mj — 1.29 and TTR = 1.00. Contours of the time-averaged streamwise 
velocity component are shown on the z — 0 plane. The plate is positioned at Xp / D — 2Q and at various radial positions, 

Vp!D. 




(c) The plate positioned at ypjD = —2. 


(d) The plate positioned at ypjD = — 1. 


Figure 9. The Md — 1.50 nozzle operates at Mj — 1.29 and TTR, = 1.00. Contours of turbulent kinetic energy are 
shown on the z — 0 plane. The plate is positioned at XpjD — 20 and at various radial positions, ypjD. 
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Figure 10. Predictions compared with measurement for the Mj — 1.29 and TTR, = 1 jet from the S1VIC016 nozzle with 
D — 0.0508 m at RjD — 100 and varying angle The plate is located at XpjD — 4 and ypjD — —2. These conditions 
correspond to F = 0.124. 




Figure 11. Predictions compared with measurement for the Mj — 1.29 and TTR = 1 jet from the SMC016 nozzle with 
D — 0.0508 m at RjD — 100 and varying angle 4^. The plate is located at XpjD — 4 and ypjD — —4. These conditions 
correspond to F = 0.027. 
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Figure 12. Predictions compared with measurement for the Mj — 1.29 and TTR, = 1 jet from the S1VIC016 nozzle with 
D — 0.0508 m at RjD — 100 and varying angle The plate is located at XpjD — 4 and ypjD — —6. These conditions 
correspond to F = 0.011. 




Figure 13. Predictions compared with measurement for the Mj — 1.29 and TTR = 1 jet from the SMC016 nozzle with 
D — 0.0508 m at RjD — 100 and varying angle 4^. The plate is located at XpjD — 4 and ypjD — —8. These conditions 
correspond to F = 0.006. 
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Figure 14. Predictions compared with measurement for the Mj — 1.29 and TTR, = 1 jet from the S1VIC016 nozzle with 
D — 0.0508 m at RjD — 100 and varying angle The plate is located at XpjD — 4 and ypjD — —10. These conditions 
correspond to F = 0.004. 




Figure 15. Predictions compared with measurement for the Mj — 1.29 and TTR = 1 jet from the SMC016 nozzle with 
D — 0.0508 m at RjD — 100 and varying angle 4^. The plate is located at XpjD — 10 and ypjD — —2. These conditions 
correspond to F = 0.309. 
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Figure 16. Predictions compared with measurement for the Mj — 1.29 and TTR, = 1 jet from the S1VIC016 nozzle with 
D — 0.0508 m at RjD — 100 and varying angle The plate is located at XpjD — 10 and ypjD — —4. These conditions 
correspond to F = 0.066. 




Figure 17. Predictions compared with measurement for the Mj — 1.29 and TTR = 1 jet from the SMC016 nozzle with 
D — 0.0508 m at RjD — 100 and varying angle 4^. The plate is located at XpjD — 10 and ypjD — —6. These conditions 
correspond to F = 0.028. 
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Figure 18. Predictions compared with measurement for the Mj — 1.29 and TTR, = 1 jet from the S1VIC016 nozzle with 
D — 0.0508 m at RjD — 100 and varying angle The plate is located at XpjD — 10 and VpjD — —10. These conditions 
correspond to F = 0.010. 






Figure 19. Predictions compared with measurement for the Mj — 1.29 and TTR = 1 jet from the SMC016 nozzle with 
D — 0.0508 m at RjD — 100 and varying angle The plate is located at XpjD — 20 and ypjD — —2. These conditions 
correspond to F = 0.618. 
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Figure 20. Predictions compared with measurement for the Mj — 1.29 and TTR, = 1 jet from the S1VIC016 nozzle with 
D — 0.0508 m at RjD — 100 and varying angle The plate is located at XpjD — 20 and ypjD — —6. These conditions 
correspond to F = 0.056. 
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Figure 21. Predictions compared with measurement for the Mj — 1.29 and TTR = 1 jet from the SMC016 nozzle with 
D — 0.0508 m at RjD — 100 and varying angle The plate is located at XpjD — 20 and ypjD — —10. These conditions 
correspond to F = 0.020. 


24 of 24 


American Institute of Aeronautics and Astronautics 


